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ABSTRACT 

The  performance  of  the  three-shell  spherical  head  model 
versus  the  performance  of  the  realistic  head  model  is  investi¬ 
gated  when  solving  the  inverse  problem  with  a  single  dipole, 
in  the  presence  of  noise.  This  is  evaluated  by  inspecting 
the  average  dipole  location  error  when  applying  a  spherical 
and  a  realistic  head  model,  for  1000  noisy  scalp  potentials, 
originating  from  the  same  test  dipole  and  having  the  same 
noise  level.  The  location  errors  are  obtained  utilizing  a  local 
linearization,  which  is  validated  with  a  Monte  Carlo  simula¬ 
tion.  For  27  electrodes,  an  EEG  epoch  of  one  time  sample 
and  spatially  white  Gaussian  noise  we  found  that  the  im¬ 
portance  of  the  realistic  head  model  over  the  spherical  head 
model  reduces  by  increasing  the  noise  level.  -Keywords: 

EEG  dipole  source  analysis;  Gaussian  noise;  realistic  head 
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The  spherical  head  model  is  constructed  as  follows. 
Through  the  27  electrode  nodes  located  at  the  scalp  sur¬ 
face,  used  in  the  finite  difference  volume  conductor  model,  a 
best-fitting  sphere  is  constructed  with  a  radius  R.  A  radial 
projection  of  the  electrodes  on  the  surface  of  the  sphere  is 
performed  to  obtain  the  electrode  coordinates  on  the  spher¬ 
ical  head  model.  The  radii  of  the  outer  shell  of  the  scalp, 
skull  and  brain  a reR,  §§R  and  f ^R,  respectively. 

When  applying  the  spherical  head  model,  the  electrode 
positions  need  to  be  obtained  to  construct  the  best-fitting 
sphere.  In  contrast,  when  applying  a  realistic  head  model, 
the  segmented  M  R-image  needs  to  be  constructed,  the  elec¬ 
trode  positions  need  to  be  obtained  and  the  numerical  pre¬ 
processing  needs  to  be  performed.  It  is  clear  that  the  latter 
is  more  demanding. 

2.2.  The  inverse  problem 

The  E EG  presented  at  the  input  of  the  inverse  solver,  with 
l  electrodes  (27)  and  for  a  single  time  instant  can  be  given 
by  Vin  e  IRlxl.  The  three  location  parameters  and  the  three 
components  of  the  dipole  are  obtained  by  finding  the  global 
minimum  of  the  residual  energy  (RE): 


RE  =  ||Vin-Vmodel||2, 

with  V modei  e  Rlxl  the  electrode  potentials  obtained  by  the 
forward  evaluation  in  the  inverse  problem  and  ||.||  the  L2 
norm. 

The  RE  indicates  the  fraction  of  energy,  which  cannot  by 
modeled  by  the  dipole.  The  Nelder-Mead  simplex  method 
is  used  to  find  the  optimal  position  ropt-  For  each  dipole 
position  the  optimal  components  dopt  are  found  in  the  least- 
squares  sense. 

2.3.  The  construction  of  AR 

For  a  given  test  dipole,  with  dipole  position  r0  and  compo¬ 
nents  d0,  the  forward  problem  is  solved  in  a  realistic  head 
model,  yielding  a  set  of  average  referenced  scalp  potentials. 
From  these  noiseless  set  of  potentials,  noisy  sets  of  potentials 
are  generated  by  adding  noise  values  to  each  of  the  poten¬ 
tials.  These  noise  values  can  bedescribed  as  a  Gaussian  with 
zero  mean.  Using  the  noiseless  forward  calculated  potentials 
in  the  realistic  head  model,  the  inverse  problem  is  solved 
in  a  spherical  head  model  yielding  dipole  position  rsphere 
and  components  dsphere-  A  local  linearization  of  the  for¬ 
ward  problem  as  presented  in  [10,11]  with  the  spherical  head 
model  at  these  dipole  parameters  is  then  performed.  This 
gives  us  the  opportunity  to  solve  for  the  deviations  of  the 
dipole  positions  Arsn^seere.  as  a  function  of  the  noise  values 
at  the  electrodes.  These  results  are  obtained  with  a  matrix 
multiplication,  which  is  faster  than  solving  the  inverse  prob¬ 
lem  iteratively.  Furthermore,  i  =  1. . .  1000  corresponds  with 
one  of  the  1000  noisy  potential  sets  with  noise  values  coming 
from  the  same  distribution  having  a  given  standard  deviation 
corresponding  with  a  given  noise  level.  The  location  error  due 
to  noise  and  due  to  the  spherical  head  model  then  becomes: 
^.sphere  =  11^*-  +  rsphere  _  ro||f  where  II  •  II  is  the 
Euclidian  norm.  For  a  1000  noisy  potentials  sets,  we  obtain 
the  average  location  error  E(Arsphere).  Next  we  calculate 
the  deviation  of  the  dipole  coordinated  A  r\eal ,  utilizing  a  re 
alistic  head  model  and  due  to  noise.  Notice  that  we  assume 
that  the  inverse  solver,  applying  the  realistic  head  model,  is 
an  unbiased  estimator.  A  local  linearization  of  the  forward 
problem  with  the  realistic  head  model  at  the  test  dipole  po¬ 
sition  r0  and  orientation  d0,  is  then  performed.  The  loca¬ 
tion  error  dueto  noise  then  becomes:  A  r\eal  =  ||Ar^e  i||, 


where  Ar^lL,  i  is  the  deviation  from  the  test  dipole  posi¬ 
tion  for  given  rioise  values.  For  the  same  1000  noisy  potential 
sets  as  in  the  spherical  head  model  casewethen  calculatethe 
average  location  error,  E{Arrieal). 

When  E(Arsphere)  is  much  larger  than  E(Arrieal),  for  a 
certain  noiselevel,  then  it  is  worthwhile  using  the  more  effort 
demanding  realistic  head  model  instead  of  the  spherical  head 
model.  On  the  other  hand,  when  E{Arlphere)  has  about  the 
same  value  as  E(Arleal),  then  it  does  not  matter  whether  a 
realistic  or  a  spherical  head  is  used.  Therefore,  thedifference 
between  the  two  average  location  errors  A  R  =  E{Arsphere)- 
E(Arleal)  is  introduced.  We  suggest  that  a  AR  smaller  than 
5  mm  indicates  that  for  a  given  dipole  and  for  a  given  noise 
level  both  models,  the  spherical  and  the  realistic  head  model 
perform  equally.  On  the  other  hand  when  AR  is  larger  than 
5  mm  for  a  given  dipole  and  for  a  given  noise  level  we  argue 
that  in  this  case  the  realistic  head  model  performs  better 
than  the  spherical  model. 

2.4.  Validation  of  the  linearization 

A  first  simulation  is  carried  out  to  evaluated  the  local 
linearization  for  one  test  dipole.  We  have  calculated 
E(Arsphere)  and  E(Arrieal)  for  a  number  of  noise  levels  and 
for  1000  noisy  potential  distributions,  utilizing  the  local  lin¬ 
earization.  On  the  other  hand  we  have  calculated,  for  the 
same  noisy  scalp  potentials,  the  inverse  problem  in  a  spher¬ 
ical  and  realistic  head  model.  This  approach  is  often  called 
a  Monte-Carlo  simulation.  The  average  location  errors  are 
derived  and  compared  with  the  ones  found  utilizing  the  local 
linearization. 

2.5.  AR  for  27  electrodes  and  Gaussian  noise 

In  thesecond  simulation,  427  test  dipole  positions,  located  in 
a  coronal  slice  containing  the  vertex  electrode  'Cz',  are  used. 

For  each  test  position  a  dipole  along  thex-,  y-,  and  z-axis  is 
generated  yielding  1281  test  dipoles.  The  x-axis  is  oriented 
from  the  left  towards  the  right  ear,  the  y-axis  is  oriented 
from  the  back  of  the  head  towards  the  nose,  and  the  z-axis 
is  oriented  from  the  bottom  of  the  head  towards  the  vertex 
electrode  'Cz'. 

Further  more  the  grand  average  over  1281  dipoles  of  A  R, 
represented  by  A(AR),  is  also  given  as  a  function  of  the 
noise  level.  Notice  that  A(  )  calculate  the  average  for  all 
test  dipoles  while  E( •)  calculates  the  average  for  all  noisy 
potential  distributions  associated  with  one  test  dipole. 

To  havea  more  quantitative  measure  for  a  given  noiselevel, 
the  number  of  test  dipoles  are  counted  which  have  a  A R 
larger  than  the  chosen  5  mm.  This  gives  us  the  opportunity 
to  investigate,  for  a  given  noise  level  and  for  a  large  number  of 
test  dipoles,  whether  it  is  still  worthwhile  utilizing  the  more 
computational  demanding  realistic  head  model,  instead  of 
the  spherical  head  model. 

3.  RESULTS 

3.1.  Validation  of  the  local  linearization 
The  results  obtained  with  the  local  linearization  are  vali¬ 
dated  utilizing  Monte  Carlo  simulations  for  one  test  dipole. 

The  solid  line  in  figure  1  presents  the  average  dipole  location 
error  E{Arrieal)  obtained  with  a  local  linearization  of  the 
forward  problem  with  thed  nic  theadh  octete,  Q^on)JCHn]2IPk]3]32K 
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noise-level 


Figure  1.  The  solid  line  (-)  and  the  dashed  line  (-  -),  depict  the 
average  dipole  location  error  as  a  function  of  the  noise  level  ob¬ 
tained  with  the  local  linearization  of  the  forward  function  when 
utilizing  a  realistic  and  spherical  head  model,  respectively.  The 
(o)  and  (x)  present  the  average  dipole  location  errors  obtained 
with  a  Monte-Carlo  simulation,  for  noise  levels  0.2,  0.14,  0.1, 
and  0.05,  utilizing  the  realistic  and  spherical  head  model,  re¬ 
spectively.  The  difference  between  the  average  location  error 
obtained  with  the  spherical  and  realistic  head  model  or  AR  for 
the  noise  level  of  0.1  is  illustrated  with  the  arrow. 


an  average  location  error.  The  'o'  depicts  the  average  loca¬ 
tion  error  obtained  from  a  Monte-Carlo  simulation  for  the 
noise  levels  0.05,  0.1,  0.14  and  0.2,  utilizing  the  realistic  head 
model.  The  'x'  depicts  the  average  location  error  obtained 
from  a  Monte-Carlo  simulation,  utilizing  a  spherical  model. 
One  can  noticethat  the  average  location  errors  obtained  with 
the  Monte-Carlo  simulation  correspond  well  with  the  values 
obtained  through  local  linearization,  and  this  even  for  larger 
noise  levels.  The  AR  =  E(Arsiphere)  -  T(Ar[eai),  for  noise 
level  0.1,  is  marked  with  the  arrow  in  figure  1.  It  can  also  be 
observed  that  by  increasing  the  noise  level  the  value  of  AR 
decreases. 

3.2.  AR  for  27  electrodes  and  Gaussian  noise 
The  first,  second,  third  and  fourth  row  in  figure  2  present 
the  value  AR  =  E{Arlphere)  -  E(Arleal)  for  a  noise  level 
equal  to  0,  0.05,  0.1  and  0.2,  respectively  as  a  function  of  the 
dipole  position.  The  dipole  orientation  is  along  the x-axis.  A 
contour  with  A  R  equal  to  5  mm  is  also  given.  For  noise  level 
zero  AR  equals  E(Arsphere).  By  increasing  the  noise  level, 
the  A  R  values  decrease  for  all  test  dipoles.  Hence  the  areas 
where  A  R  is  larger  than  5  mm  decrease  when  increasing  the 
noise  level. 

The  grand  average  over  all  the  test  dipole  positions  of  A  R 
is  presented  in  the  column  marked  with  A(AR)  in  table  1.  It 
can  be  observed  that  by  increasing  the  noise  level  the  grand 
average  of  A  R,  decreases  from  5.5  mm  to  1.8  mm. 

I  n  table  1,  the  percentage  of  A  R  values  larger  than  5  mm 
are  also  given,  for  the  applied  noise  levels.  Notice  that  60% 
of  the  test  dipole  have  a  A  R  larger  than  5  mm  for  the  noise¬ 
less  case.  Evoked  potentials  and  epileptic  spikes  coming  from 
the  same  focal  source  are  often  averaged,  to  reduce  the  back¬ 
ground  EEG.  For  the  averaged  spikes  the  noise  level  is  typi¬ 
cal  situated  around  0.1.  For  this  noise  level  13.4%  of  the  test 
dipoles  have  a  A  R  larger  than  5  mm.  This  number  decreases 


Figure  2.  The  first,  second,  third  and  fourth  row  illustrate  A  R  = 
E(Arsphere)  —  E(Ar\eal)  for  noise  levels  equal  to  0,  0.05,  0.1 
and  0.2,  respectively.  A  contour  with  AR  equal  to  5  mm  is 
also  given.  The  configuration  with  27  electrodes  and  Gaussian 
noise  is  used.  The  results  for  dipoles  oriented  along  the  x-axis 
are  depicted. 
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Table  1.  The  grand  average  over  the  1281  test  dipole  for  AR,\s 
presented.  Further  more,  the  relative  number  of  test  dipoles 
yielding  a  A R  larger  than  5  mm,  is  also  given.  The  noise  levels 
0,  0.05,  0.1  and  0.2  are  considered.  27  electrodes,  and  Gaussian 
noise  is  applied. 


noise  level 

SNR 

A(AR)(mm) 

(|(A R  >  5  mm)(  %  ) 

0.0 

oo 

5.5 

“60U 

0.05 

20 

4.0 

25.5 

0.1 

10 

3.1 

13.4 

0.2 

5 

1.8 

7.6 

to  7.6%  for  a  noise  level  of  0.2,  which  corresponds  with  an 
unaverage  epileptic  spike.  The  higher  the  noise  level,  theless 
important  becomes  the  usage  of  the  realistic  head  model  over 
the  spherical  head  model. 

4.  CONCLUSION 

The  aim  of  this  paper  was  to  evaluate  the  performance  of  the 
three-shell  spherical  head  model  versus  the  performance  of 
the  realistic  head  model,  in  solving  the  inverse  problem  with 
a  single  dipole  in  the  presence  of  noise.  In  our  simulations, 
the  average  location  errors  E(Ar*phere)  and  E(Arrieal)  are 
obtained  through  a  local  linearization.  As  control  we  found 
for  one  dipole  that  the  values  of  E{Arlphere)  and  E(Arffal) 
obtained  through  local  linearization  are  in  good  agreement 
with  the  ones  obtained  when  utilizing  a  Monte  Carlo  ap¬ 
proach.  When  utilizing  27  electrodes,  an  epoch  of  onetime 
sample  and  Gaussian  noise,  it  was  found  that  by  increasing 
the  noise  level  from  0  to  0.2,  that  the  grand  average  over 
all  the  test  dipoles  of  AR  decreased  from  5.5mm  to  1.8mm, 
respectively.  For  a  noise  level  of  0.1,  corresponding  with  an 
averaged  epileptic  spike,  13.4%  of  the  test  dipoles  have  a  A  R 
larger  than  5mm.  For  a  noise  level  of  0.2,  corresponding 
with  an  unaveraged  epileptic  spike,  less  than  10%  of  the  test 
dipoles  have  a  AR  larger  than  5mm.  The  importance  of  the 
realistic  head  model  over  the  spherical  head  model  reduces 
by  increasing  the  noise  level. 
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